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Abstract 

The low lying spectrum of the 0(n) effective field theory is calculated in the delta- 
regime in 3 and 4 space-time dimensions using lattice regularization to NNL order. 
It allows, in particular, to determine, using numerical simulations in different spatial 
volumes, the pion decay constant F in QCD with 2 flavours or the spin stiffness p 
for an antiferromagnet in d = 2 -|- 1 dimensions. 



1 Introduction 



The low energy phenomena in systems with spontaneously broken symmetry are governed 
by the dynamics of the Goldstone bosons. This can be described by an effective field theory, 
and the calculations could be performed by chiral perturbation theory (ChPT) [1] |2]. 
The effective action contains low energy constants (LEG) determined by the underlying 
microscopic theory. The physical quantities can be systematically expanded in powers of 
momenta, or (as in the case of our interest here) inverse box size. In numerical simulations 
one can place the system into a space-time box of size Lt x and study the dependence 
of different quantities on the box size. Comparing the data with theoretical predictions 
one can determine the corresponding LEC's. 

There are important cases when the order parameter of the spontaneous symmetry 
breaking is an 0(n) vector. These include QCD with two light quarks (4d 0(4) case), and 
antiferromagnetic layers (3d 0(3) case). 

The different regimes of such systems in a finite box have been systematized by Leutwyler 
[3]. In particular, for the case of no explicit symmetry breaking (zero quark mass in QCD) 
one can distinguish the e-regime ("cubic geometry"), where Lt ~ L^, and the 5- regime 
("cylindrical geometry") where Lt 3> l/m{Ls) ~ F'^L^^ 3> Lg. (Note that the expansion 
parameter of ChPT in this case is 1/(F^L^) and must be small enough.) 

In this work we calculate the low lying spectrum in a finite spatial box in the 0(n) 
effective theory (i.e. we consider the ^-regime) with no explicit symmetry breaking, in 3 
and 4 space time dimensions, using lattice regularization. It has been shown in [3] that 
in the leading order the spectrum is given by the quantum mechanical 0(n) rotator with 
moment of inertia O = (n — 1)/2 ■ F^L^ (the "angular momentum" being the 0{n) isospin). 
In the case of QCD the constant F is the pion decay constant. The next-to-leading (NL) 
order term of the expansion in 1/(F^L^) has been calculated in [4j. In the calculation one 
considered an infinitely long lattice in the time direction, separating the spatially constant 
slow modes and the fast modes, and integrating out the latter. The resulting effective 
Lagrangian for the slow modes is then that of an 0(n) rotator, with a modified moment 
of inertia. The NLO correction turned out to be large: at Lg = 2.5 fm it is still ~ 40%. 
Therefore it was important to calculate the NNLO term. This has been done recently for 
the 4d, 0(4) case by P. Hasenfratz [5j by a method similar to the one used in [1], except 
that in [5] dimensional regularization (DR) was used. The NNLO term contains two new 
LEC's, Ai and A2, and is estimated to be —5% at Ls = 2.5 fm. 

Since the calculation using DR with the time- dependent slow modes was quite involved, 
we have chosen to calculate the same quantity by a different method, following the cal- 
culation of the small- volume mass gap in 2d 0{n) non-linear sigma-model by Liischer, 
Weisz, Wolff using lattice regularization, and we considered both the 3d and 4d cases 
for general 0(n). 
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As mentioned above, the effective theory has also been successfully applied in condensed 
matter physics. In particular one can perform numerical simulations in the microscopic 
theory of the spin | antiferromagnetic Heisenberg model and measure different quantities 
(like staggered susceptibility, etc.) with an impressive precision [71 18]. Comparing these 
with the results of the effective field theory one can obtain the LEC's (spin stiffness, etc) 
to high accuracy of order 10~^ (for a recent paper see [9]). Given the high accuracy, in this 
case even a small NNLO term could have an important effect. 

The rest of the paper is organized as follows. Section [2] discusses the different terms 
in the effective action needed to this order, section |3] recapitulates the method of [B] to 
obtain the perturbative expansion for the mass gap. In section 4 we describe the numerous 
checks of our calculations, and section 5 gives the general expression for the mass gap. The 
numerical values and the renormalization of the couplings for d = 4 and d = 3 dimensions 
are given in sections 6 and 7, respectively. Section 8 contains our conclusions. Some further 
details of the calculations are delegated to the Appendices. 



2 The effective action 



We consider the following effective action in 3 and 4 space-time dimensions 

A = A2 + A4 + ..., (2.1) 

where the leading term is 

A2 = 7^y^d,S,-d,S,, x = (a:o,x), 8^ = 1. (2.2) 

Here denotes the standard forward lattice derivative, and the A4 part containing 4- 
derivative terms will be specified later. Note that we do not consider here an explicit 0(n) 
symmetry breaking term. 

The non-linear sigma model given by A2 is non-renormalizable in c? > 2 dimensions 
hence one has to consider A as an effective low energy action, with corresponding (infinitely 
many) low energy constants (LEG). The mass gap (and the energy of lowest states) in a 
box can be expanded in inverse powers of the box size Lg, and to a given order in 
l/Lg only finite number of LEG's appear. For our case of NNLO corrections we need only 
terms up to four derivatives in d = 4. 

The action (12. ip in d = 3 describes e.g. the low energy behaviour of the antiferromag- 
netic spin I quantum Heisenberg model, while in d = 4, and n = 4 it is the effective action 
of QGD with two massless quarks. 
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The dimensionless lattice coupling Aq in eq. fl2.2p is expressed through the dimensionful 
lattice bare parameters of the corresponding theories as 



= — , for = 3 , (2.3) 

poa 

A2 = -i^, ford = 4, (2.4) 

where po and Fq denote the bare spin stiffness and bare pion decay constant, respectively, 
and a is the lattice spacingQ We parametrize the spin vector with the "pion" fields 



AoTf, ^1 - A^Tf^J , 7f= (7r\7r2,...,7r"-i). (2.5) 

The pion fields reflect the perturbative fluctuations around the magnetization axis. We 
expand in the vf-fields keeping the volume Va = Lf x finite assuming a cylindrical 
geometry, Lf ^ Lg. As long as the volume is finite there is a slowly moving global 
mode which corresponds to the direction of the magnetization. This mode is treated 
nonperturbatively in the the path integral [TO] . 

The partition function takes the form 

Z = Xll f difr,6(yJ2'^J e-^=«['^l (2.6) 

where 

Acff [n] =A[tt]+ A^casure [vf] + A.ero [vf] . (2.7) 

The first term is the action of (12.11) (expressed in terms of vf- fields). The second term is 
caused by the measure from the change to vf-variables 



^measurc[vf] = lu (l - A^TT^) ^ . (2.8) 

X 

The irrelevant factor Af and the last term together with the delta function arise when we 
separate the global zero mode and integrate it out in the path integral, 

Aero[v?] = -{n - 1) In 5^ (1 - A^vf^^ (2.9) 

X 

As a consequence of the delta function in (12. 6p we must leave out the p = (0, 0) zero mode 
in the momentum decomposition of vf. 



^In most expressions we use the convention "a = 1", but restore the a factors when it is useful for 
understanding. 
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For the d = 4 case at NNLO we also need the A4 term in eq. fl2.ip . The most general 
4-derivative interaction consistent with the lattice symmetries is given by (see e.g. [TT] ) 



where 



and 



A 
A 
A 

A, 
A 



i=2 



X 
X 



2 

X J J 



2 

XJ } 



Y^Y,{d,S,.d,S,f--j^(A'-^^+2AP] , 



X fj, 



d + 2 



X n 
4(56) _ 4(1) 



A\ ^ — ^ ^ d^dij^Sx ■ dydySx- 



(2.10) 

(2.11) 
(2.12) 
(2.13) 

(2.14) 
(2.15) 

(2.16) 

(2.17) 
(2.18) 



We use the standard forward [d) and backward {d*) lattice derivatives. Some comments 
are in order here: 



• The interaction A\ is redundant at the order of our calculation. Transforming 
variables in the path integral as S — )■ (S + anS)/|S + anS| the leading action exhibits 
a contribution proportional to ^4^''. Choosing the parameter a we can absorb the 
contributions of this interaction [T2] . 

• In ChPT of QCD with Nj = 2 flavours (in Minkowski space) the standard notation 
for the low energy constants is k p3]. Following [12] for the 0{n) case (in Euclidean 
space) we use the convention g^^ /4: for the couplings. They are related by 

9? = -^lu 9T = -^l2. (2.19) 

The minus sign comes from going from Minkowski to Euclidean space. Note also 
that the convention for numbering the operators is different. 
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• The operators 4 and 5 are absent in dimensional regularization. They are needed 
to restore Lorentz symmetry on the lattice and their coefficients are fixed by this 
requirement. The subtracted pieces stem from the fact that we construct first the 
symmetric, traceless 4-index tensors, and take the sum X^^^mmmm [H]- 

• According to equation f l2.10p we subtract a term proportional to the leading action 
A2 from each of the 4-derivative interactions. The coefficients c^*-' serve to remove 
the power-like divergence from the contribution of the corresponding operator. 
The subtracted operators renormalize then multiplicatively. This is discussed later. 



3 Correlation function 



We extract the mass gap from the correlation function 

= y^^(^^o,^-Syo,y)^ t=\yo-Xo\, (3.1) 

where Vg = L'jr^ is the spatial volume. We follow closely the method of ref.[n] developed 
to obtain the small-volume mass gap m{Ls) for the 2d 0{n) nonlinear sigma-model. 

We apply perturbation theory at fixed finite cylindrical volume Lt x and following 
[6j we impose free boundary conditions in the time direction^. This guarantees that we 
project onto states with zero isospin (and zero total momentum) at the two boundaries. 
Except for the ground state all these states ("scattering states") have energies > Att/Ls. 
The correlation function drops off as 



C{t) = Ae-""^^^^' + O (e-^J*) , (3.2) 



where m{Ls) = Ei — Eq is the mass gap, the difference between the lowest energies in the 
/ = 1 and / = isospin sectors. For small Aq the mass gap is m{Ls) oc Ag/K hence for 
Aq ^ Vs/Lt the system is nearly completely polarized - the spins fluctuate only slightly 
around the direction of the total magnetization. Therefore for flxed Ls,Lt one can use 
perturbation theory in Aq to calculate the correlation functioijfl 



C{t) = l + J2{^IY^'CA^) . (3.3) 



i=0 



L 



Inserting the expansion of the mass gap m{Ls) into eq. (13. 2 p one concludes that the coef- 
flcients Ci{t/Ls) are (up to exponentially decreasing terms) polynomials of order i + 1 in 



^with periodic b.c. in the spatial directions 

•^The notations are analogous to those of ref. [6^, except that there a symmetric set-up with xq — —t, 



yo = T, t = 2t is used. 
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the time t 



(3.4) 



As seen from eq. (13. 2p (neglecting the exponentially small terms in t) the expansion of 
logC(t) is linear in t, the higher order terms in t should cancel. To obtain the mass gap 
to NNLO only the coefficients Coo, Coi, Cio, Cn, and C21 are needed, while Cu, C22, and 
C23 are useful to check the absence of t"^ and terms in logC(t). 

Observe that the use of free b.c. in the time direction is essential here. Using periodic 
b.c. also in this direction, the transitions /—)•/ + 1 between states of higher isospins with 
energies Ei ^ l{l + n — 2)Xl/{2Vs) would also contribute to the perturbative expansion of 
C{t), making the method unpractical for the determination of the mass gap. 

Note that the method of ref. [6] is essentially an e-regime expansion in a very elongated 
cylindrical volume, since the spins are strongly correlated over the whole length Lt. The 
correlation length C,t{Ls) = 1/iti{Ls), defined for an infinitely long cylindrical volume be- 
comes much larger for Aq — )■ than any finite Lt . For the "truly 5-regime" calculation 
one should study the dynamics of the spatially constant slow modes, which is described 
by the quantum mechanical rotator [3llll|5|. Of course, the two approaches should lead to 
the same result for the mass gap, but one meets different technical difficulties in these two 
approaches. 



4 Procedure and checks 

It is straightforward to derive the Feynman rules and to write down the corresponding 
Feynman diagrams. After that we separate the different contributions analytically from 
each graph. This step is practically the same for d = 2 and d > 2. Evaluating the diagrams 
numerically provides a good check for separating the different powers of t. There are also 
other checks for the final results: 

• The consistency relations mentioned in section [3] are satisfied. 

• The NLO result is known for general dimensions [1]. 

• The NNLO result for d = 2 can be compared to the results of Liischer and Weisz 

[ni[i5] 

• We also considered the same problem with Dirichlet-free boundary conditions. (This 
should obviously give the same result for the mass, but the corresponding transition 
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amplitudes are different.) We checked the result numerically for = 2, including the 
4-derivative contributions. 

• Since the 4-derivative contributions were not calculated previously on the lattice for 
the 2d case, we compared the corresponding contributions to direct Monte-Carlo 
simulations for c? = 2, = 3 at sufficiently small Aq values. 

• We solved the problem for n = 2 parametrizing the fields as S(x) = (cos sin 
including the 4-derivative contributions. 



5 The mass gap at finite lattice spacing 



Restoring the lattice spacing a the mass gap reads 



2Lei 



i=2 



+ 0{Xl) (5.1) 



The n-dependence of the coefficients is given by 

C2 = C21 + 022^1 , 

C3 = C31 + 032^1 + 033^1 , 

4) = rfg+45ni, 2 = 2,3,4,5. 



where we introduced the abbreviation 



Til = n — 1 . 



(5.2) 
(5.3) 
(5.4) 

(5.5) 



The above coefficients are expressed as lattice sums over spatial momenta and are given in 
appendix [Bl 



6 Results for d = A 



The expressions occurring in the previous section depend on the ratio a/Lg only. In four 
space-time dimensions one has 

C2k{a/Ls) = C2k0 + C2klj2 + ^ [jlj ^ ^ = 1,2, (6.1) 

C3fc(a/Lj = C3fco + C3fcl-^ + C3fc2-^ + C3fc3-^log-^ + C ( -^J , /c = 1, 2, 3 , (6.2) 

42(a/L.) = Dg + E«^ + (9(^] , ^ = 2,3,4,5; A: = l,2. (6.3) 
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In eq. flS.ip these coefficients are multiplied by Ag = I/FqO^ and Ag = I/FqC^, respectively. 
Therefore the omitted terms correspond to lattice artifacts. The resulting terms with 
power-like and logarithmic singularities in a will be absorbed by the renormalization of the 
lattice bare parameters. The numerical values of the coefficients are listed in Tables [T]l3] in 
appendix O 



6.1 Discussion on the numerical values of appendix [C 



We have calculated the lattice sums in appendix [B] to a quadruple precision for different 
lattice sizes and ffited the a/Ls-dependence. In order to find out a reliable set of terms 
one can vary the minimum value of Lg and include fewer or more higher powers like a^/L^, 
a^/L^, etc. The coefficients should be stable under such variations, and we can estimate the 
precision of the values, which was 10 digits or better. Note, that we know the coefficients of 
the logarithmic terms from DR, hence fixing these parameters in the ffi allows to determine 
the remaining ones to better accuracy. 

The coefficient C211 = — C221 is the shape coefficient for a three dimensional cubic 
box, given in [12]. The NLO finite size effect was ffist calculated in [1]. 



The logarithms originate from the double sums in flB.6p and f lB.7p . To obtain the 
leading Lg-dependence we expand them for small ki, k2, 

1 «2 

'''^2V^,^ (51 + ^2 + s,)s,s,s, + • • • ^^-^^ 
K1K2 



C32 



1 // 



where we expanded in Si = sinh/i,. (See appendix [X] for notations.) The leading singular- 
ities of the double sums are related to those in the sum 

X 

_ 1 ^„ (1 + e-^i-'^2-M3)(cosh/ii - 1) 1^ 
~ 4:V^ ^ (1 - e-'^i^^a-^s) sinh/iisinh/i2sinh/i3 4 ^ 8^ ^ 

— Y" d + .... (6.6) 

^^s (S1 + S2 + S3)S1S2S3 



Therefore, one has 



C31 = 2A + . . . , (6.7) 
C32 = -2A + . . . . (6.8) 
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According to the calculation in DR p] the logarithmic part of A is given by 

J dx DoGix)Gixr = -y ll^^ + • • • ' (6-9) 

where q = 0.837536910696 (cf. (IC.ip ). The coefficient of the logarithmic term in C31 is, 
therefore, —0.0353584296400 in agreement with the direct fits given in Table [2] in appendix 
O 



6.2 Renormalization for = 4 

The additive renormalization of the 4-derivative operators in eq. fl2.10p are removed by 
setting 

c^^^ = (^S + ^i^S) , ^ = 2,3,4,5. (6.10) 

The singular part of the a-dependence of the mass gap can be removed by the renormal- 
ization of the bare lattice parameters as 



1 


1 


^2 




(i) 
94 


(i) 

= 9 m 



(6.11) 



gS log aM = gS log aMi, z = 2,3, (6. 12) 

where M is a scal^. Like in DR we introduce the individual scales M2 and M3. Note that 
g^^ and gf '' do not need renormalization to this order. The coefficients in eq. (16.111) are 
expressed through the infinite volume limits of coefficients in eqs. (15.21) . (15.31) as 

bi = -C2(0) = G(l, 0) - niG(O) = - (n - 2)G(0) , (6.13) 
62 = -C3(0) + 2c2(0)2. (6.14) 

Their numerical values for d = 4 are 

61 = 0.029933390231 - 0.154933390231 ni , (6.15) 

62 = 0.001425585601 - 0.019893440100 ni + 0.024004355409 nj . (6.16) 

Note that there are several nontrivial relations that should be satisfied to be able to ab- 
sorb the singular cutoff dependence into the renormalized couplings - these are satisfied 
numerically to the expected accuracy, providing an additional check. 

The renormalization of g^^^ and gf^ agrees with the result of ChPT in dimensional 
regularization [2], taking into account eq. (I2.19p ^. On the lattice (with Euclidean action) 



*In order to avoid confusion with DR wc do not use A or /i. 

^In ref. only the n — A result is given but it is straightforward to restore the general n dependence. 
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we have 



(2) 
^41 



(3) 

9ai 



1 /14 

2ni 

167r2 V 3 

1 8 

^16^3 ■ 



(6.17) 
(6.18) 



Note that, in contrast to the DR, the pion decay constant Fq renormahzes on the lattice. 
The NLO coefficient bi in (16. lip has been calculated earlier in [1]. 

ChPT in the p-regime with lattice regularization was studied earlier by Shushpanov 
and Smilga [16] for the 4d 0(4) case, who obtain, besides other 1- and 2- loop results, 
also the renormalization of Fq to NLO. However, their result disagrees with ours - they 
obtain (for n = 4) bi = —2G{0), i.e. the —l/{2d) term of eq. (16.131) is missing. It is 
easy to verify its presence for the 0(2) case in the same way as done in [16] (from the 
current-current correlator) when one uses the parametrization by the angular variable. In 
this case the action is given by ^(1 — cos{\o^f^(f){x))) while the conserved current is 

A^{x) = A|^^ sin(Ao9^0(a;)). In this case one gets bi = —l/(2d), as expected. 



6.3 The moment of inertia 



The isospin dependence of the lowest excitations up to (and including) NNLO corrections 
are given by the rotator spectrum [TTl [5], i.e. by 

^i = ^^^^, / = 0,1,2,.... (6.19) 

where G is the moment of inertia. We have also calculated the E2 — E0 gap with the method 
presented here, and found agreement with this expectation. The moment of inertia is given 
by 



e = L^F^ 



1 + —^7:^0.225784959441 {n - 2) 
F^Lj 

^ ' -0.0692984943 + 0.0101978424 n) 



-^0.007071685925 [(3n - 10) log M2L, + 2n log M3L,] ^q_20) 

(4) 

[-0.55835794046 (n + 1)] 



F^LV 



(5) 

J^[0.55771822866 - 1.11639602502 + C . 
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The scales M2 and M3 are related to the corresponding scales Ai, A2 in dimensional regular- 
ization (cf. eq. f l2.19p ). However, to get this relation, and the values of the coefficients gf* 
and g^p (needed to restore Lorentz symmetry) one needs to relate the lattice regularization 
to the dimensional one. This step remains still to be done. The resulting uncertainty is a 
term const/ (F^L^). The remaining terms are in agreement with the result of ref. [5] where 
the calculation was performed in DR for the n = 4 case. 



7 Results for d = ?> 



In three space-time dimensions the coefficients of (I5.2p - (l5.4p depend on a/L^ as 

C2fc(a/Ls) = C2fco + C2fci^ + . . . , /c = l,2, (7.1) 

(X 

C3fc(a//^s) = csfco + csfci-^ + C3fc2-^ + • • • , A; = 1,2, 3, (7.2) 



42(a/L,)=D2 + 0(^^J, ^ = 2,3,4,5; A; = 1,2. (7.3) 
The corresponding values are given in Tables |1J16] in appendix O 



a 



7.1 Renormalization and the moment of inertia for d = 3 

The 4-derivative interactions only contribute to order l/L^. (Note that the corresponding 
couplings g^l^ are dimensionful, in contrast to c? = 4). They are also not needed for the 
renormalization here. Because the theory is non-renormalizable, it is expected that they 
are necessary to absorb divergences at higher orders. The spin stiffness renormalizes as 

l.l(l + ^L + Ji, + ...). ,7,4) 

The corresponding coefficients are still given by eqs. fl6.13p . fl6.14p while their numerical 
values are 

61 = 0.0860643431920 - 0.252731009859 ni , (7.5) 

62 = 0.0102138509611 - 0.0659002864141 m + 0.0638729633447 (7.6) 



and the moment of inertia is 

e = pL: 



n — 2 n — 2 /I 
0.310373220693 —0.000430499941 + O 



pLs P^Ll \Ll 



(7.7) 
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The coefficient of the NNLO term has been estimated a long time ago from finite 
temperature simulations in the spin | quantum Heisenberg model [S] and the result cited 
there is not consistent with our calculations. The discrepancy of parameters obtained 
from those finite-temperature data with other high precision measurements (about two per 
mille, but statistically significant) has been observed already in ^ and it was attributed 
to unaccounted finite-temperature corrections. A new measurement [IB] of the staggered 
susceptibility at much lower temperatures agrees with our very small coefficient of the 
oc term. 



8 Summary 

We calculated the mass gap in the 0{n) effective field theory in a cubic spatial box of 
size Ls for 3 and 4 space-time dimensions to NNLO, using lattice regularization. The 
renormalization of the bare lattice couplings is performed, however, the connection of the 
4-derivative couplings with the MS scheme is not established yet. This affects only the 
const/F^L^ term in the NNLO term in d = 4 (relevant to QCD with two massless quarks), 
not the logarithmic terms, where we reproduced the result of p]. Note that the effect of 
including a small explicit symmetry breaking (a small quark mass) has been done to LO 
in [3], and to NLO recently in [12]. In d = 3 (relevant for 2-1-1 dimensional spin ^ quantum 
Heisenberg model) our calculation is complete to NNL order since the 4-derivative operators 
do not contribute to this order. 
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A The propagator 

We expand the leading action in the vr-fields and write it as 

X x,y 

The kernel p{x, y) decomposes into the spatial and time direction 

d-l 

pix,y) = ^ps{Xf,-yf,) + po{xo,yo)- (A.2) 
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The spatial kernel Ps{x) is up to the sign the standard one- dimensional lattice Laplace 
operator 

Ps(x) = 25a;,o - - 5x,L,-i, X = 0,1, . . . ,Ls - I, (A.3) 

and periodically continued for x > Ls and x < 0. It is convenient to write pQ{xo,yo) for 
free boundary conditions in matrix form 



Po 



-1 


2 


-1 .. 


. 











.. 


. -1 


2 -1 








.. 


. 


-1 1 



Consider the partial Fourier transform of the kernel 

p{xo, yo, k) = ^ p(xo, Vo, X - y)e'^(''-y^ = pq{xo, yo) + S^g^y^PsCk) 



(A.4) 



(A.5) 



where 



d-i 



p,(k) = 5^4sin2(A;,/2) = 4sinh2(/x(k)/2). 



i=l 



(A.6) 



This equation defines /i(k) which is denoted by p in the following. It will be useful to 
separate the k = contribution in the propagator 



G{x, y) = G{xq, yo, X - y) = ^ G(xo, yo] pi^^W^^"" ^'^ 



Y9oixo,yo) + Gi{x,y) = Go{xo,yo) + Gi{x,y). 



(A.7) 



The 1-dimensional massless propagator is the inverse of po in flA.4p with the zero mode left 
out (the boundaries are at t = ±T) 



go{xo,yo) = -2^^o - yo\ + 



, 2^ , TiT+1] 



and 



2(2T + 1)' " 3(2T + 1)' 



(A.9) 



The primed sum symbol means that the k = zero mode is excluded. The function 
G{t,t'] p) is the massive propagator for the Id case with free b.c: 



G{t,t'-p) 



2sinh/i(l -e-2/^(2T+i) 

{ 



^1-*-*')}. (A.IO) 
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This function is exponentially small for |t — t'| ^ L^. In the limit T — )■ oo and when both 
time arguments are finite the propagator is simplified, 

p-Ml*-*'l 

G{t,t'-fi)^-——. (A.ll) 
2 smh n 



The oc dependence of a Feynman graph for C{t) comes from the time-like part 
Go{xo,yo) in dMl). 



B Analytic results 

B.l Contributions from the leading action 

We introduce the notations 



The coefficients read: 



k 
k 



k 



C3i{a/Ls 



1 



1 



2^2 12 



5i 



12K 



-105i - 3^3 + 6^2 + 12^0 



+ ^^i(C3 + 2£i-25i)-^ 



+ 



1 — V // 

^ kik2 ^ 



6-^1-^2-/^3) sinh. sinh sinh /i3 



[< 
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C32(a/L,) = ^(52 - ^3) + ^(15^2 + 25i + 10^0 



1 



1 ^ — \ // 

^ (T 



24(rf - 1 



-(Ci + 65i-6C3)Ci 



96K^ ^ (1 _ e-^i-^'2~''3) sinh fii sinh /i2 sinh 1x3 

kik2 



_4g-2Ati-2At2 _j_ ^g-2^ii-2^2-2Ai3j (B 7) 

(The double prime in the sum over ki, k2 means the condition that ki 7^ 0, k2 7^ 0, and 

kg = ki - k2 ^ 0.) 



C33(a/L^) 



1 



1 1 

2 sinh /i 



(B.8) 



B.2 Contributions from the next-to-leading action 

In this section we give the contributions of the 4-derivative interactions, without the con- 
tribution of the subtracted terms proportional to c*^*-* in equation fl2.10p . The effect of these 
terms is discussed in the renormalization chapter. The contributions of the interactions 
2,3,4 depend only on the expression 



^0 — d^d^Gy 



1 1^,1- e-^' 
Vs Vs ^ sinh /i 



(B.9) 



They are given by 



7(2) 

^31 

j(3) 
^31 

7(4) 
^31 



-2O0 

2d 
^ d + 2 



(B.IO) 
(B.ll) 
(B.12) 



"32 

(3) 
32 



d 



-no 



7(4) 

*32 



d + 2 



d 



(B.13) 
(B.14) 

(B.15) 



The contribution of the fifth operator is more complicated. Using the labels introduced in 
section |2] we write 



d 



(5) 
3i 



d + 2 



24f ) + S"'^ 



^3j 



i = 1,2. 



(B.16) 
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We obtain 

~ Vs (l + e-M)3 +^^Z. 2sinhV^^ coskO,(B.17) 



A5a) _ 
"32 



1 1 "^,(3 + 6 ^)e 1 v^' cosh/i ^ 2 /-DION 



1 1 



4? = --, (B.20) 



j(5c) _ 
"31 



1 1 „ , (3 - Se-'' - 4c-2m + e-^'')e-'' 
"t7 t7 2^ 



(1 + 6-/^)2 



d-1 



(l + e-M)2smh/. 2^(l-cosk0 



«=i 



1 1 



ci-l 



V; ^ 2 sinh^ 



^e*'(l-coski) 



— (cosh/x — 1)' 



, (B.21) 



i=l 



j(5c) _ 
"32 



2(1 + e-/^)^ 



d-l 



1 ' (1 - 2e-^ - e-2^)e-'^ . 
~ K ^ n i^-M^2«^^H„ 2^^^ ~ coski)^ 



(1 + e~^)^ sinh 



1 / cosh /i 



K 4;^ 2 sinh'' 



d-l 



^e^^'(l -coski) 



i=l 



— (cosh fjL — ly 



(B.22) 
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C Numerical values 



C.l Numerical values for d = 4 

The values for the coefficients occurring in equation (16. ip and the following are collected 
in Tables [T112] where we introduced the notations 



q = 0.837536910696 , r = 0.9764866840 , 
Vq = -0.03163189123 , = 0.01546809528 . 



k 


C21fc 


C22fc 





-0.029933390231 


0.154933390231 


1 


0.225784959441 


-0.225784959441 



Table 1: Coefficients at NLO, d = A. 



k 


Csifc 









0.000366430100 


0.001342713582 


0.024004355408 


1 


-0.013517018599 


0.083480277057 


-0.069963258459 


2 


0.0961483532 


-0.051187951 


0.108634522 


3 


-0.035358429639 


0.035358429639. 






Table 2: Coefficients at NNLO, = 4. 



k 


^(2) 


^3k 


^(3) 
^3k 


^3k 


^(4) 
^3k 


^(4) 
^3k 


^(5) 
^3k 




1 


1 
2 


-2q 


5 
4 


-Q 





4q 
3 


vo 




2 


-1 





1 
4 







2q 
3 




3 1 



Table 3: Coefficients for the 4-derivative contributions at NNLO, d = 4. 



C.2 Numerical values for d = 3 

The values for the coefficients occurring in equation (17. ip and the following are given in 
Tables ma 
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k 


C21fe 


C22fe 





-0.086064343192 


0.252731009859 


1 


0.310373220693 


-0.310373220693 



Table 4: Coefficients at NLO, d = 3. 



k 




C32fe 


Cask 





0.004600291377 


-0.021104227057 


0.063872963344 


1 


-0.053424134767 


0.210306009764 


-0.156881874998 


2 


0.095901036182 


-0.192232572305 


0.096331536123 



Table 5: Coefficients at NNLO, d = 3. 



k 


^(2) 
^3k 


^(3) 


^(4) 
^3k 


^(5) 
^3k 


1 


2 
3 


4 
3 





0.038725856392 


2 


-1 


1 
3 





0.008715670880 



Table 6: Coefficients for the 4-derivative contributions at NNLO, d — 3. 
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